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Abstract 

A narrow dip structure has been observed at 1.9 GeV/c^ in a study of diffractive 
photoproduction of the 37r"'"37r~ final state performed by the Fermilab experiment 
E687. 



1. Introduction. 



The Fermilab Experiment 687 collaboration has collected a large sample 
of high-energy photoproduction events, recorded with the E687 spectrometer 
[1][2] during the 1990/91 fixed-target runs at the Wideband Photon beam- 
line at Fermilab. Although the experiment is focussed on charm physics, a 
very large sample of diffractivcly photoproduced light-meson events was also 
recorded. This paper reports on a study of the diffractive photoproduction 
of the 37r'^37r~ final state and the observation of a narrow dip in the mass 
spectrum at 1.9 GeV/c^. 

2. E687 spectrometer 

In E687, a forward multiparticle spectrometer is used to measure the inter- 
actions of high-energy photons on a 4-cm-thick Be target. It is a large-aperture, 
fixed-target spectrometer with excellent vertexing, particle identification, and 
reconstruction capabilities for photons and 7r°'s. The photon beam is derived 
from the Brcmsstrahlung of secondary electrons of ~ 300 GeV endpoint en- 
ergy, which were produced by the 800 GeV/c Tevatron proton beam. The 
charged particles emerging from the target are tracked by a system of twelve 
planes of silicon microstrip detectors arranged in three views. These provide 
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Fig. 1. Distribution of 37r"''37r~ invariant mass after applying a cut on the total 
energy deposited in the calorimeters with respect to the total energy in the 
spectrometer. 

high-resolution separation of primary (production) and secondary (charm de- 
cay or interaction) vertices. The momentum of a charged particle is determined 
by measuring its deflections in two analysis magnets of opposite polarity with 
five stations of multiwire proportional chambers. Three multicell threshold 
Cerenkov counters are used to discriminate between pions, kaons, and pro- 
tons. Photons and neutral pions are reconstructed by electromagnetic (EM) 
calorimetry. Hadron calorimetry and muon detectors provide triggering and 
additional particle identification. 



3. Event Selection 



Pions are produced in photon interactions in the Be target. The data acqui- 
sition trigger requires a minimum energy deposition in the hadron calorimeters 
located behind the electromagnetic calorimeters and at least three charged 
tracks outside the pair region. The microstrip system and the forward spec- 
trometer measure the Gvr final state (in this paper, 6n refers to the Stt+Stt" 
state) with a mass resolution a = 10 MeV/ c^ at a total invariant mass of about 
2 GeV/c^. It is required that a single six-prong vertex be reconstructed in the 
target region by the microstrip detector, with a good confidence level. Such 
a requirement rejects background due to secondary interactions in the target. 
Exclusive final states are selected by also requiring that the same number 
of tracks be reconstructed in the magnetic spectrometer. The six microstrip 
tracks and the six spectrometer tracks are required to be linked, with no am- 
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Fig. 2. Tranverse momentum squared distribution showing the diffractive peak. 



biguity in associating the microstrip and spectrometer tracks. Events with 
particles identified by the Cerenkov system as definite electrons, kaons, or 
protons, or as kaon/proton ambiguous are eliminated and at least four out of 
six particles have to be positively identified as tt^. Particle identification is 
tested by assuming that one or two out of the six tracks is a K^, by comput- 
ing all two-track invariant mass combinations, and verifying that there is no 
evidence of a peak at the K* or at the mass. We ehminated final states with 
7r°'s by rejecting events with visible energy in the electromagnetic calorime- 
ters that was not associated with the charged tracks. A cut in this variable 
{Eca,i/EQ„ < 5%) is applied on the calorimetric neutral energy normalized 
to the six-pion energy measured in the spectrometer. The distribution of the 
six-pion invariant mass after these cuts is shown in Fig. 1. The plot shows a 
structiuc at 1.9 GeV/c^. In the following, only the 67r mass region around this 
structure will be analyzed. 



For diffractive reactions at our energies, the square of the four-momentum 
transfer t can be approximated by the square of the total transverse momen- 
tum P| of the diffractively produced hadronic final state. Using this defini- 
tion, the distribution of diffractive events. Fig. 2, is well described by two 
exponentials: a coherent contribution with a slope be = 54 ±2 (GeV/c)^^ con- 
sistent with the Be form factor [3] and an incoherent contribution with a slope 
hi = 5.10 ±0.25 (GeV/c)-2. Taking only events with P| < 0.040 GeVVc^, we 
evaluated a contamination of about 50% from nondiffractive events. This in- 
coherent contribution shows no structure in the 1.2—3.0 GcV/c^ mass range. 
Fig. 3. The diffractive mass distribution was obtained by subtracting this con- 
tribution, parametrized by a polynomial fit, and dividing the yield by the 
detection efficiency. 
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The detection efficiency was computed by modeling diffractive photopro- 
duction of a mass M, using tlie experimentally found slope b^, and simu- 
lating the decay 6tt according to phase space [4]. There is no thresh- 
old or discontinuity for the efficiency, particularly in the region of the dip 
structure. At 1.9 GeV/c^, the (self-normalized relative) efficiency A varied as 
dA/A/dM^T, = 0.15/ GeV/c^. The efficiency and the efficiency- corrected dis- 
tribution of the six-pion invariant mass for diffractive events, in the mass range 
1.4—2.4 GeV/c^, arc shown in Fig. 4. There no evidence, albeit with large com- 
binatorial backgrounds, for resonance substructure, e.g., p° — > tt'^tt" , in the Gtt 
data below Mq^^ — 2.0 GeV/c^, either at the mass region of the dip or in nearby 
sidebands. Similarly, the efficiency or acceptance exhibited no threshold, edge, 
or discontinuity over the entire mass region observed, when the six pion state 
was simulated as a sequence of decays of intermediate two-body resonances, 
for example + ai {p^t^^) + {p^t^~) ~^ (7r"'"7r~7r"'") -|- (7r"'"7r~7r~), even under 
extreme assumptions of full longitudinal or transverse polarizations for the 
initial state. 

The presence of a dip at 1.9 GeV/c^ was verified by several checks of 
the systematics. Diffractive photoproduction of D^D^ pairs or the associated 
production of plus a charm baryon at low t (where the decay products of 
the other charm meson or charm baryon are missed) followed by the 6 vr decay 
of the D° or are estimated by Monte Carlo simulation to be negligible 

contributions. It was also checked that demanding more stringent cuts (i.e., 
requiring that all six particles be identified as tt^, applying a sharper cut 
on the calorimeter neutral energy, or subtracting the incoherent contribution 
bin by bin) increases the statistical errors without significantly affecting the 
behavior shown in Fig. 3. 

A three-parameter polynomial fit was performed, solid line in Fig. 4, to ex- 
plore the hypothesis that any structure in this distribution is a statistical fluc- 
tuation. The normalized residual distribution, evaluated for each 10-MeV/c^ 
bin, is good in the full invariant mass range (Fig.5: left), with the exception 
of the interval centered at 1.9 GeV/c^, the region of the claimed dip, where a 
poor ~ 10~^ confidence level interval (Fig. 5: right), is obtained, making it 
highly unhkely that the observed dip is a statistical fluctuation. Incoherently 
adding a Breit-Wigner to the flt does not improve the flt quahty much, dashed 
lines in Fig. 4 and Fig. 5. 



4. Fitting the six-pion invariant mass distribution 



Because of the narrow width, the E687 spectrometer mass resolution, a ~ 
10 MeV/c' at 2 GeV/c^, was unfolded by applying the method described in 
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Fig. 3. Distribution of Svr^Svr^ invariant mass in the 1.2—3.0 GeV/c^ mass range: 
coherent plus incoherent contribution. Dotted distribution: incoherent 
contribution. 
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Fig. 4. Acceptance-corrected distribution of 37r"'"37r~ invariant mass for diffractive 
events after subtracting incoherent contribution. Sohd hne: second-order 
polynomial fit. Dashed line: polynomial fit with incoherently-added Breit-Wigner. 
Upper dot line: relative detection efficiency (arbitrary normalization). 

Ref. [5]: the experimentally observed data distribution r{x) and the unfolded 
mass distribution a(a;) are related by a(a;) ~ r{x) — 0.5(j'^-r{x) , where r(x) 
is the second derivative with respect to Mqtt for the observed distribution. This 
relationship results from applying a Fourier transform and approximating the 
resolution function by g{x) = exp{— ^j^^ ), which is marginally different from 
a Gaussian. A fit similar to the following one is used to obtain r{x)" from the 
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Fig. 5. Normalized residual (left) and confidence level (right) distributions. 

Solid line: second-order polynomial fit.Dashed line: Breit-Wigner plus second-order 

polynomial fit. 

unfolded data. The data after unfolding are shown in Fig. 6. 

The dip structure at 1.9 GeV/c^ has been characterized by a two-component 
fit, adding coherently a relativistic Breit-Wigner resonance to a diffractive 
continuum contribution. The continuum probability distribution Fjs{M) has 
been modeled after a Jacob-Slansky diffractive parameterization [6], plus a 
constant term cq 

-13 
p M—Mq 



The Jacob-Slansky model is based on probabilities, rather than amplitude 
phenomenology, and specifically has no prescription as to how the complex 
phase of the continuum distribution varies with mass. The Jacob-Slansky am- 
plitude fjs{M) is assumed to be the purely real {(pjs = 0) square root of the 
probability function Fjs{M). For the fit, a relative phase factor e^'^, indepen- 
dent of mass, and a normalizing factor multiplied a relativistic Breit-Wigner 
resonance term, giving the overall amplitude 

— M Fp**^ 

A{M) ^ fjs{M) + a 



M2 - M2 + iM^r 



Fit results are shown in Table 1 and in Fig. 6 for a fitted mass range from 
1.65 to 2.15 GeV/c^, symmetric with respect to the dip. Quantities shown 
are the mass and width of the resonance, the amplitude ratio a,./ f]s{Mj.) 
between the Breit-Wigner function and the Jacob-Slansky continuum, the 
relative phase and the x^/dof. The Jacob-Slansky parameters for that one 
mass range are also given. Fit values show consistent evidence for a narrow 
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Table 1 

Fit results for a mass range from 1.65 to 2.15 GeV/c? 



M., (GeV/c2) 


1.911 ±0.004 


r (MeV/c^) 


29± 11 


ar//js(Mr) 


0.31 ±0.07 


(j) (deg.) 


62± 12 


xVdof 


1.1 


Mo (GcV/c^) 


1.49 ±0.02 


Co (GcV/c2)-i 


0±50 


ci(GeV/c2)i-" 


960 ± 80 


13 (GeV/c2) 


0.5 ±0.3 


a 


1.8 ±0.2 



resonance at M,. = 1.911 ± 0.004 ± 0.001 GeV/c^ with a width T = 29 ± 11 ± 
4 MeV/c^, where the errors quoted are statistic and systematic, respectively. 
The fit values shown in Fig. 6 and represented by the parameters of Table 1 
are stable with acceptable x^/dof over varying mass ranges from 1.65 to 2.3 
GeV/c^. We quote as systematic error the sample variance of the fit values due 
to our choice of fit mass range. The quality of the fit deteriorates somewhat 
as the upper limit of the fit for this simple model is extended from 2.1 to 
3.0 GeV/c^. However, the only fit parameter that is affected is the width, 
which varies from 29 ± 11 MeV/c to 40 ± 20 MeV/c. 



The dip structure is reminiscent of what was observed in e+e~ annihilation 
[7] [8]. The mechanism by which a narrow resonance may interfere destruc- 
tively with a continuum, or a broad resonance, could be similar to the one 
described in [9], in a different context. Vector qq hybrids are predicted at ~ 
1.9 GeV/c^ according to the fiux tube model [10] [11]. A hybrid is expected 
to have a small, but not vanishing, e"''e~ width and a total width constrained 
by decay selection rules [12]. Vector glueballs are expected at higher masses, 
according to the bag model and to lattice calculations [13] [14]. Narrow reso- 
nances have also been predicted near the A^'A'" region, but have never been 
conclusively found [15] [16] [17] [18] [19]. A A^A^ threshold effect might also pro- 
duce a downward step in the amplitude, followed by a recovery [20] [21]. These 
resonances should also be observed as a steep variation in the nucleon time-like 
form factor [22] [23]. 
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Fig. 6. Acceptance-corrected distribution of Stt^Stt" invariant mass for diffractive 
events. The mass resolution has been unfolded. Fit parameters are listed in 
Table 1. 



5. Conclusions 



The diffractive photoproduction of 3n'^3n~ has been studied by E687. 
Evidence has been found for a narrow structure near MgTr = 1-9 GeV/c^. If 
this dip is characterized as the destructive interference of a resonance with the 
continuum background, then the parameters of this resonance would be 
= 1.911 ± 0.004 ± 0.001 GeV/c2, with T = 29 ± 11 ± 4 MeV/c^. Such a 
resonance could be assigned the photon quantum numbers (J^'-' = 1 ) and 
G=+l, 1=1 due to the final state multiphcity. There is httle understanding of 
the specific mechanism responsible for this destructive interference. In order 
to facilitate additional phenomenological studies, the data points of Fig. 6 are 
available [24]. 
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